Microscopic Calculation of Spin Torques in Disordered Ferromagnets 
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Effects of conduction electrons on magnetization dynamics, represented by spin torques, are cal- 
culated microscopically in the first order in spatial gradient and time derivative of magnetization. 
Special attention is paid to the so-called /3-term and the Gilbert damping, a, in the presence of 
electrons' spin-relaxation processes, which are modeled by quenched magnetic (and spin-orbit) im- 
purities. The obtained results such as a ^ (3 hold for localized as well as itinerant ferromagnetism. 

PACS numbers: 72.25.Ba, 72.15.Gd, 72.25.Rb 



Recent activity in spintronics to manipulate nanoscale 
magnetization, such as magnetization reversal, by electric 
current started with the concept of spin-transfer effect 
introduced by Slonczewski pj and Berger When a 
current flows in a background magnetization n varying 
in space, a slight directional mismatch (with n) arises 
in the electron spin polarization such that spin angular 
momentum of a spin-polarized current (spin current) is 
transferred to the magnetization. This effect is expressed 
as a torque on n, called spin-transfer torque; for slowly- 
varying n, it takes the form, — (t> s -V)n, where «s is a 
velocity characterizing the spin-transfer rate 0, 0,1MB 0- 

Existence of another type of current-induced spin 
torque, of the form, — (3n x (v s ■ V)n, was noted re- 
cently H, 0, 0]. This torque is perpendicular to the 
spin-transfer torque, and is called nonadiabatic torque 
or the /3-term. Zhang and Li showed on phenomenologi- 
cal grounds that it arises as a result of spin relaxation of 
conduction electrons . Barnes and Maekawa suggested 
its close connection to the Gilbert damping, a, and pro- 
posed a relation a = f3 |l(J . While the importance of the 
/3-term in magnetization dynamics is now well-recognized 
1^.13. HfUllj . its microscopic derivation is a current issue. 

Similar torques are known also in multilayer [l2l Il.'i| , 
domain wall 0, 0, 0] , and other 0, systems. For 
example, for a rigid domain wall, a torque which acts 
as a force and may thus be called momentum-transfer 
torque, was noted |lj,[l6|]- This is a spatially oscillating 
torque due to electron reflection, found by Waintal and 
Viret 0], which has the same alg ebraic form as the /3- 
term but is spatially nonlocal [lj|. This torque can be 
important for magnetic configurations varying rapidly in 
space whereas the /3-term is relevant to slowly-varying 
configurations. We focus on the /3-term in this Letter. 

Very recently, a microscopic study of this subjest was 
undertaken by Tserkovnyak, Brataas and Bauer (TBB) 



based on the Boltzmann equation 20] . They showed that 
a = f3 for an itinerant 'single-band' ferromagnet, and 
emphasized its peculiar magnetization dynamics. Their 
analysis is, however, still phenomenological as to the 
spin-relaxation term which was introduced by hand. 

In this Letter, we present a microscopic calculation 
of spin torques, especially the /3-term and the Gilbert 
damping, in a localized (s-d) as well as an itinerant 
(Stoner) models for ferromagnetism. We introduce spin- 
relaxation processes into the (conducting) electron sys- 
tem by magnetic (and spin-orbit) impurities causing spin- 
flip and spin-dependent scatterings. Our result shows 
that both transverse and longitudinal relaxation pro- 
cesses contribute to a and /3, and that a =^ /3 in general. 

We first consider a localized (s-d) model. It consists of 
localized d-spins, S, and conducting s-electrons, which 
are coupled via the s-d exchange interaction 



H sd = -M 



J dr n(r)-a(r). 



(1) 



Here we put S = Sn with a unit vector n pointing in 
the direction of spin |2l|, <x(r) = cJ(r)crc(r) represents 
(twice) the s-electron spin density, with c) = (c|,c{) be- 
ing electron creation operators, er the Pauli spin-matrix 
vector, and M = J S &S with J s d being the s-d exchange 
coupling constant. The total Hamiltonian of the system 
is given by H tot = H s + H c \ + H sd , where H s and H e i 
are for localized d-spins and s-electrons, respectively. 

The dynamics of magnetization, — n [2l|. will be de- 
scribed by the Landau-Lifshitz-Gilbert (LLG) equation 



n = 70-Hefl x n + aoh x n + t' e 



(2) 



where "foH c g and «o are an effective field and a Gilbert 
damping constant, respectively, both coming from Hs 
and not from processes involving s-electrons. Effects of 
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conducting s-electrons are contained in the spin torque 

toi(r) = -g-^(r) = Mn(r) x (a(r)) n . e ., (3) 

which comes from H s d- (Here a 3 is the volume per local- 
ized spin.) The calculation of spin torque is thus equiv- 
alent to that of s-electron spin polarization, (<7(r)) n .c, 
or precisely, its perpendicular projection (<x_i_('"))n.e. to 
n |22j, in such nonequilibrium states with current flow 
and spatially varying magnetization (for the /3-term), or 
with time-dependent magnetization (for Gilbert damp- 
ing). Here and hereafter, (• • • ) n . e . represents statistical 
average in such nonequilibrium states. 
Generally, the spin torque is expressed as 

t e i = aon + (a- V) n + bo (n x n) + n x (6 • V) n (4) 

in the first order in time derivative and spatial gradients. 
The a-term includes the spin-transfer torque, — (u s -V)n, 
where u , = — (a 3 /2eS)j s (e > being the elementary 
charge) |21| is a quantity which may be called as the 
spin-transfer velocity. (The ao-term just renormalizcs 
the magnitude of spin, S — > S — (a 3 /h)ao on the left- 
hand side of the LLG equation, or in front of the kinetic 
("Berry phase") term in the spin Lagrangian 7].) The 
6-term is the central issue in the present study. We also 
focus on 60 1 which gives damping of the Gilbert type, 
with dimensionless damping constant a = — (a 3 /hS]bn. 
The b-term is expressed as — (3n x (u s -V)n (in t' el ) |2lJ 
in the literature |9j, [l(J, |2£j , which defines the constant (3. 

The spin torque (@J corresponds to the s-electron spin 
polarization (<xj_) n . e . \22\ given by 

[boh + (b- V) n — ao (n X h) — n X (a -V) n] . (5) 

To calculate the coefficients, a M and 6 M , we follow TBB 
[2o| (see also Ref.[23|) and consider a small transverse 
fluctuation, u = (u x ,u y ,0), \u\ <C 1, around a uniformly 
magnetized state, n = z, such that n = z + u. In the 
'unperturbed' state, n — z, the s-electrons are described 
by the Hamiltonian 

Wo = J2( £k ~ £F ^ C ka C ka + V imp, (6) 

A; ci- 
and have a spin polarization (<r)o = psZ. Here p s = 
n^—ni, n„ = fc Fo ./67r 2 , efo- = h 2 kp a /2m = sp+aM, and 
V[ mp is the impurity potential specified later (Ea.l|I5fr). 
The subscript a =T,| corresponds, respectively, to a — 
+1, —1 in the formula (and to a =|, t or —1, +1) . In the 
presence of u(r,t) = u(q, uj) e l (i' r ^ UJt ) ^ the s-electrons 
feel a perturbation (note that H cl + H sd = Hq + Hi) 

Ut = -M^(^(r Cfc .T.((i lW )e- fa * 1 (7) 

ha 

and acquires a transverse component [22j 

(^x(^))n.e. = M [ dt' f dr' X f(x-x')u^(x')(8) 

J — oo J 



in their spin polarization, in the first order in u. Here, 
x = (r,t), x' = (r',t'), and \^ is the transverse spin 
susceptibility in a uniformly magnetized state, (a, (3 = 
x, y specify the transverse componentsC and summing 
over (3 is implied.) Writing in Fourier components, 

(af(q,u)) n . e . = Mx a /(q,^ + iO)u^(q,uj), (9) 

we expand x± (<7, w + *0) as 

X f («,a; + tO)-xf (0,0) 
= -j^[i(b-q ~ b ui)S a [3 + i(a-q - a u!)e a p} (10) 

up to the first order in q and to. (Here e a p is an antisym- 
metric tensor in 2D, with e xy = 1.) Below we will see that 
Mx± (O,0) = p s S a (3- Therefore, (<x)„. . = p s n+ (a-±) n . c ., 
with {<T±_) n .c. given by 

[b u + (b-V)u - a (z x it) - z x (a- V)it] , (11) 

which coincides with (jSJ) within the present accuracy (i.e., 
neglecting u x u = 0(u 2 )), and leads to the spin torque 
(J3J ■ The problem thus reduces to the calculation of the 
coefficients, a M and 6 M , in Ea. l|TUj l 2J- 

We consider a 3D electron system with = h 2 k 2 /2m, 
which is affected by the impurity potential represented by 

M mp = uY.Hr-R^+u^Srcrdir-R^in) 

i 3 

in the first-quantization form. The first term describes 
potential scattering, and the second describes spin scat- 
tering by impurity spins, Sj. We take a quenched average 
for the impurity spin direction as 

Sfsf = M«/9x(S {a .^^ X ^ (13) 
I Sj (a,/3 = z) 

as well as for the impurity positions, Ri and R'j , as usual. 
The electron damping rate is then given by 

7(T = -L = Txn iU 2 v a + Trn s u 2 (2Sl v s + ~S%v a ) (14) 

in the first Born approximation, where n\ (n s ) is the con- 
centration of normal (magnetic) impurities, and v„ — 
mkYa/2iT 2 Ti 2 is the density of states (DOS) at e^a- Wc 
can include the spin-orbit impurities, of scattering am- 
plitude m s . .(fe x k')-cr, by adding n s , u 2 (k x fe')| to 
n s u 2 S 2 , in all expressions in the present study. 

We assume that 7 CT <C £f<t and 7 CT <M, and calculate 
a and /3 in the lowest non-trivial order in x a = "fa/^e-po 
and 7<r/M, both being collectively denoted as 7. We will 
see that a and /? are quantities of 0(7). For a perturba- 
tive estimate in 7, only the first-order correction is suffi- 
cient to the a a - (and a^-) vertex, since the corresponding 
'diffusion' ladder has a large mass ~ M. 
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Gilbert damping : We first study the w-linear terms 
in the uniform (q = 0) part of the transverse spin suscep- 
tibility, X±(Q = 0,cj + iO). Contributions up to 0(7), 
represented in Fig. 1 , are given, in Matsubara form, by 

X*f(0,iujx) = -T^^idap + iaea^ifa, (15) 



where 



'■Po 



(16) 



with Xc = Xa(i£ n ;iuJ\) = Y.k G ks{i£ n + iuJ\)G ka {ie n ), 
Gka(z) = (z - Ek + £Fa + ^erSgn^mz)) -1 . The vertex 
correction is associated with a factor T = riiu 2 — n s u 2 S 2 , 
where the S z scattering changes sign whereas S x and S y 
scatterings cancel each other. After analytic continuation 
(iuj\ — > uj+iO) and extracting the w-linear terms, we have 



bo = 



M 2 h 
2tt 



5]Re{^(0;0)-^ 2 )(0;0)}, (17) 



at absolute zero, T = 0. Here the upper labels, (1) and 
(2), on ip a (and x<y below) indicate the analytic continu- 
ations, G(ie n + iuj\)G(ie n ) — » G R G R and G R G A , respec- 
tively 25]. 

By an explicit evaluation of the fc-integrals, we have 
xl 1} (0;0) = (»7r/2M)[i/_ + iy - i{v_^/M)] + 0( 7 2 ), 
xS 2) (0; 0) = -{i-K/2M)[au+ - iy + i{u +1± /M)] + 0( 7 2 ), 
and thus b = -%(v + ^ + — v-^-)/A+(-K%to/S) (v\ -v 2 _). 
Here v± = if ±1/4., 7 ± = ( 7 t±7i)/2, and y = v^-v^x^. 
Using Ea. H14|) . we finally obtain 



2S 2 y^ i +S\(v 2 + v 2 ) 



S 



(18) 



As expected, only the spin (and spin-orbit) scattering 
contributes to a, and the potential scattering (~ n;u 2 ) 
does not thanks to the cancellation between selfenergy 
and vertex correction |26| . 

By a similar analysis, we obtain ag — —hp s /2. 

For iuj\ — 0, ip a is independent of a. Thus, from 
Ea. lflT)) . x" /3 (0, 0) is proportional to 5 a f3- Explicit calcu- 
lation shows that Xj^O'O) = (Ps/M) $af3 as used above. 

/3-term : We next examine the q-linear terms in the 
presence of current flow under static n(r) = z + tt(r). 
We produce a current-carrying state by applying a d.c. 
electric field E, and calculate a linear response of <r" to 





FIG. 1: Transverse spin susceptibility up to 0(7). The dotted 
line with a cross represents scattering by impurities, either 
non-magnetic or magnetic (including spin-orbit). The thick 
(thin) solid line carries Matsubara frequency ie n +icu\ (ie„)- 



E. In a similar way to the Kubo formula for electrical 
conductivity j2^| , one can derive 



(^(g))„.e. - x?{q)Ei, 



(19) 



= (20 ) 



dre l ^ T (T T a a (q,T) J,). (21) 



Here J — — e J2k VkC ka c ka i s the total current (vi~ = 
hk/m), and the average (• • • ) is taken in the thermal 
equilibrium state determined by Hq + Ti-i with lu = (in 
Eq.(JZJ|). To have a non-vanishing contribution, we need 
to extract u 13 and qj both in first order, and put 

K?(q,iu x ) = -eMKg{iw x )q j u f, (q,Q). (22) 

Diagrammatic expressions for K°j^ are given in Fig. 2. 

The calculation of b (and a) as a linear response to 
E is similar to the calculation of transport coefficients, 
but requires a closer analysis to collect all the terms one 
order higher in 7. The result, up to 0(7°), is 

Eh h . 



efiE 
2m 

efiE 
2mM 



^ an a T a ^ s - 7rf ^ on a T a v s (24) 



A [(5| + Sl)v + j s + (SI Sj)u.j c ] , (25) 

where j s = a s E = jij — j± is the spin current, and j c — 
a c E = j| + j[ is the charge current, with {a c ,a s } = 
(e 2 /m)(n|r|±n|T4) being charge and spin conductivities. 
Similarly to the case of a, only the spin (and spin-orbit) 
scattering contributes to b 26]. We also see that b is 
mainly determined by j s , but also depends on j c for a 
general case of S" 2 ^ S 2 , and that 



(3 



M 



(sl + si) v+ + -(sl~si) v - 




v, 



,(26) 






FIG. 2: The E t -, qj- and it^-linear coefficient, K°f , of the 
transverse spin polarization, (a™ (q))a.e., in the presence of 
current flow (Ei) and magnetization texture (qju 13 ). The ve- 
locity vertices with Vi and Vj are associated with Ei and qj, 
respectively. Other graphical meanings are the same as Fig.l. 
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where Pj = u s /a c is the polarization of the current. 

In terms of longitudinal (t l ) and transverse (r T ) spin- 
relaxation times, obtained as 1 = 4:Trn s u 2 Sj_v+/?i and 
t t x = 27m s u 2 (Sl + S^)i^ + /ti we have 



P = 



AS 

h 

2M 



Pi)- 



P z — 



p-i 



(27) 
(28) 



where P v = v-/v + is the DOS asymmetry. We see that 
both longitudinal and transverse spin-relaxation pro- 
cesses contribute to a and (3 in general (i.e., if P v ^ 0), in 
contrast to the demonstration of Ref.[20|. For "isotropic" 
impurities with Sj_ = S 2 and thus r L = r T = r s , we have 



a 3 is + % fi 



(29) 



Thepresent results are obtained from those of Zhang and 
Li y| if we identify their r s f and no with our r s (or t t 
when P v = 0) and Mv+, respectively. 

Finally, we consider an itinerant single-band model 
such as the Stoner model treated in a mean-field ap- 
proximation. Suppose we apply a small transverse field, 
h = h qu e^ q ' r ~ UJt \ on a uniform ferromagnetic state. 
The corresponding magnetization dynamics is described 
by u a (q,uj) = — [x*± PA ( ( l> w ) ] Q ^q w i within linear re- 
sponse. Here u = (<T^_)n.o. now, and x± PA = (x± — ^) _1 
is the RPA (transverse) susceptibility, with x± given by 
Eq. (|10|l and J being the ferromagnetic exchange coupling 
constant. We write (x I ± PA )~ 1 u = —h, and expand as 

Xl 1 = Xl^o ~ XlfiX±,iX±fi U P t0 °(<li u )' where x±,n 
represents terms n-th order in q and oj. Using the self- 
consistent equation, M = Jp s , at equilibrium, we obtain 

p~ 2 [b u + (b-V)u- z x (a u + (o-V)«)] = ft.. (30) 

Within 0(u), this is consistent with the LLG equation 

h = {a 1 - V)n + n x [b' h + (b' -V)n - p 2 s h'}. (31) 

Primed quantities here are those divided by — ao = hp s /2. 
This amounts to replacing 2S in the localized model by 
p s a 3 . For "isotropic" impurities, we have a — hv + /2p s T s 
and P = h/2Mr s , and thus a ^ p even in an itinerant 
single-band model, in contrast to the result of Ref . |20j . 

We have presented a model calculation of spin torques, 
especially the Gilbert damping and the so-called /3-term, 
on the basis of two types of microscopic models for ferro- 
magnetism (localized and itinerant) and a controlled ap- 
proximation. Magnetic (and spin-orbit) impurities have 
been used as a model for spin relaxation in conducting 
electron systems. In actual systems, however, the origin 
of spin relaxation will be various depending on specific 
systems. We hope the present formalism will be useful for 
subsequent studies treating such more realistic models. 
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